Non-classical Rotational Inertia in a Two-dimensional Bosonic Solid Containing Grain 
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We study the occurrence of non-classical rotational inertia (NCRI) arising from superfluidity 
along grain boundaries in a two-dimensional bosonic system. We make use of a standard mapping 
between the zero-temperature properties of this system and the statistical mechanics of interacting 
vortex lines in the mixed phase of a type-II superconductor. In the mapping, the liquid phase of 
the vortex system corresponds to the superfluid bosonic phase. We consider numerically obtained 
polycrystalline configurations of the vortex lines in which the microcrystals are separated by liquid- 
like grain boundary regions which widen as the vortex system temperature increases. The NCRI of 
the corresponding zero-temperature bosonic systems can then be numerically evaluated by solving 
the equations of superfluid hydrodynamics in the channels near the grain boundaries. We find that 
the NCRI increases very abruptly as the liquid regions in the vortex system (equivalently, superfluid 
regions in the bosonic system) form a connected, system-spannig structure with one or more closed 
loops. The implications of these results for experimentally observed supersolid phenomena are 
discussed. 



I. INTRODUCTION 

The observation 1 ^ of non-classical rotational inertia 
(NCRI) in torsional oscillation experiments on solid 4 He 
created a great deal of interest^ in the possibility of occur- 
rence of a "supersolid" phase in which crystalline order 
and superfluidity coexist. Considerable subsequent work 
has sustained this interest. The occurrence of NCRI in 
solid 4 He has been confirmed in many experiments^— and 
the dependence of the measured NCRI fraction (NCRIF) 
on various factors, such as the method of sample growth, 
its annealing, frequency and amplitude of the torsional 
oscillator and the amount of 3 He impurities present, have 
been studied in great detail. However, despite these ex- 
tensive studies, the microscopic origin of the observed 
NCRI signal remains quite controversial. In particular, it 
is not clear whether the superfluid component in super- 
solid 4 He is distributed uniformly throughout the sam- 
ple, or present only near structural defects such as dis- 
locations and grain boundaries. A suggestion that the 
observed phenomena were related to an old theoretical 
proposal^ of superfluidity arising from mobile zero-point 
vacancies in the crystal has been contradicted by results 
of quantum Monte Carlo simulation s 11 ' 12 that show that 
the concentration of vacancies is too low to account for 
the measured NCRIF. On the other hand, the observed 
dependence^ of the magnitude of the NCRIF on the qual- 
ity of the sample (samples with higher degree of crys- 
talline order exhibit smaller NCRIF) argues in favor of a 
mechanism of superfluidity in which defects play a major 
role. An important role for the defects is also indicated 
by the observation*^ of a close correspondence between 
the onset of NCRI and an increase of the shear modulus 
of the solid, and by an enhancement of both these effects 



when the concentration of 3 He impurities is increased. It 
has been suggested^ 3 - that both the occurrence of NCRI 
and the simultaneous increase in the shear modulus arise 
from a stiffening of a network of dislocation lines, and 
that this stiffening is assisted by 3 He impurities which 
pin the dislocation lines by binding to them. 

The question of whether the experimentally observed 
supersolid behavior in 4 He can arise from the presence of 
structural defects has been investigated, on the theoreti- 
cal side, in quantum Monte Carlo studies 1 ^ - — that show 
that superfluidity can indeed occur along grain bound- 
aries and in the cores of dislocations. There have been a 
few attempts^— to calculate the macroscopic properties 
(such as the NCRIF) of a system in which superfluidity 
occurs along a network of grain boundaries or disloca- 
tion lines. Since these structural defects form disordered 
complex networks, a calculation of the rotational inertia 
of a superfluid confined in a network of irregular-shaped 
channels is necessary for understanding whether a defect- 
based mechanism can provide a consistent explanation 
of the observed results. To our knowledge, no such cal- 
culation for realistic defect network structures currently 
exists in the literature. 

In this paper we study the possibility of occurrence of 
"supersolid" behavior, similar to that observed in solid 
4 He, arising from superfluidity along grain boundaries 
in a two-dimensional bosonic system. We make use of 
a well-known mapping 22 . - — between the statistical me- 
chanics of a system of interacting vortex lines in the 
mixed phase of type-II superconductors and the zero- 
temperature quantum mechanics of a two-dimensional 
system of bosons. Specifically we consider a highly 
anisotropic layered superconductor in the presence of a 
magnetic field perpendicular to the layers. The system 
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then consists of a collection of vortex lines oriented, on 
the average, normal to the layers. The statistical me- 
chanics of these vortex lines can be mapped^r— onto 
the quantum mechanics of a two-dimensional system of 
interacting bosons at zero temperature. The mapping 
can also be used if a small concentration of columnar 
pinning centers normal to the layers is present in the vor- 
tex system. In that case, the equivalent bosonic system 
contains a small concentration of randomly located point 
pinning centers that produce a random external poten- 
tial for the bosons. In this mapping, the liquid phase of 
the vortex system corresponds to the superfluid phase of 
the bosons. The equilibrium properties of the vortex sys- 
tem, both with and without columnar pinning, have been 
investigated in a large number of theoretical^—, experi- 
mental^— and numerical^ - — studies. These results es- 
tablish, as we will show below, that the zero-temperature, 
two-dimensional, interacting bosonic system with ran- 
dom pinning centers exhibits a polycrystalline state with 
superfluidity along grain boundaries over a suitable range 
of system parameters. 

Our earlier studies^ - — of the vortex system in the 
presence of columnar pins provide us with several real- 
istic configurations of the network of grain boundaries. 
We find that some of these polycrystalline states survive 
as metastable states when the random pinning potential 
is turned off, as explained below. This case maps then 
onto the impurity-free bosonic system, which is similar 
to the experimentally studied 4 He case. We then study 
the NCRIF arising from the superfluid regions along the 
grain boundaries by numerically solving the equations of 
superfluid hydrodynamic o 33 ' 34 in the geometry specified 
by the superfluid channels in the sample. This allows us 
to determine the NCRIF as a function of system param- 
eters. At low temperatures (meaning the temperature 
of the vortex system, not that of the two-dimensional 
bosonic system), the liquid regions (superfluid regions in 
the equivalent two-dimensional system of bosons) near 
the grain boundaries are completely absent or small, and 
the NCRIF is vanishingly small. As the temperature is 
increased towards the melting temperature of the vortex 
lattice in the absence of any pinning, a kind of "pre- 
melting" occurs near the grain boundaries, so that the 
area covered by the liquid regions increases. This in- 
creases the connectivity of the network of liquid channels 
and causes the NCRIF to increase as the channels open 
up throughout the entire sample. However, the NCRIF 
remains vanishingly small as long as the liquid regions 
remain isolated from one another. Our main result is 
that the NCRIF exhibits a sharp jump to a measurable 
value of a few percent when the growing liquid regions 
percolate across the system to form one or more closed 
channels of size comparable to the size of the system. 
The behavior of the NCRIF as a function of the temper- 
ature of the vortex system is qualitatively similar to that 
seen in the experiments on 4 He. Our work, thus, shows 
that supersolid behavior similar to that observed in solid 
4 He can occur from superfluidity in grain boundaries in 



a two-dimensional system of bosons for an appropriate 
choice of parameters. 

The rest of the paper is organized as follows. In sec- 
tion UH we provide the details of the vortex-to-boson 
mapping used in our work, define the model we consider 
and describe the numerical methods used in our calcula- 
tions. The results of our study are described in detail in 
section Hill Section ITVl contains a discussion of the impli- 
cations of our results in the context of current research 
on supersolidity in 4 He. 



II. MAPPING, MODELS AND METHODS 

In this section, we first describe in detail the mapping 
between the statistical mechanics of a collection of vortex 
lines and the zero-temperature quantum mechanics of a 
two-dimensional system of bosons. We then discuss how 
existing results for the vortex system can be used to infer 
the occurrence of a "supersolid" phase with superfluid- 
ity along grain boundaries in the bosonic system under 
suitable conditions. In our study, information about the 
behavior of the bosonic system is obtained, through the 
vortex-boson mapping, from calculations carried out for 
the vortex system. At the end of this section, we provide 
some details of the model used in our study of the vortex 
system and the method of calculation we have used. 



A. Mappings and Models 

It is well established^— that the partition function of 
a system of interacting vortex lines can be written as that 
of a two-dimensional system of interacting bosons. We 
review here only the necessary details. Let us consider 
a system of N vortex lines in a type-II superconductor 
in the mixed phase, with the magnetic field in the z- 
direction. We also assume that N p columnar pinning 
centers oriented along the z direction may be present in 
the sample (this random pinning potential will be turned 
off in the numerical work described below). Let the two- 
dimensional vector Vj(z) denote the transverse position 
of the j-th vortex line in the xy-p\a,ne at z (0 < z < L 
where L is the thickness of the sample in the z-direction) . 
The free energy of the system of vortex lines has the form 
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Here, e is the tilt modulus of the vortex lines, V(r) is 
the interaction potential between two vortex lines sep- 
arated by transverse distance r, and Vd(r) is a pinning 
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potential produced by the columnar pinning centers (for 
columnar pins oriented in the z-direction, the pinning 
potential does not depend on z). We assume periodic 
boundary condition in the z-dircction, which implies that 



the positions {ij(0)} and {rj(L)} at the two ends of the 
sample must match modulo a permutation. The parti- 
tion function of the system of vortex lines at temperature 
T is then given by 
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(2.2) 



where the functional integrals are over all vortex-line configurations that satisfy the boundary conditions Tj (L) = 
r P(j) (0) f° r a ^ i> P representing a permutation of the indices 1, 2, . . . , N. 

To see the connection of this problem with the quantum mechanics of bosons, let us consider a two-dimensional 
system of N identical bosonic particles of mass to, with pairwise interactions given by the potential V(r), in the 
presence of an external impurity potential Vd{v). In the path integral representation, the partition function of this 
system at temperature T b can be written as 
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where Tj(t) now represents the position of the j-th boson at "imaginary time" r and A = H/T b . A comparison of this 
expression with that in Eq. (|2.2I) shows that the two are the same when one makes the identifications 



e ^ to, T ->■ 7i, L ->■ A = h/T b . 



(2.4) 



Thus, the thermodynamic limit, L — > oo, in the vortex 
system corresponds to the zero temperature (T b = 0) 
limit in the boson system. The tilt modulus of the vor- 
tex lines plays the role of the mass of the bosons and the 
temperature T of the vortex system plays the role of H. 
The pairwise interaction V(r) in the system of bosons is 
the same as that between two vortex lines and the exter- 
nal potential Vd(r) is also the same in the two systems. 
The interaction between two straight vortex lines parallel 
to each other is repulsive and its dependence of the trans- 
verse separation is given by V(r) oc Kq(t/X) where Kq is 
a Bessel function and A the in-plane London penetration 
depth of the superconductor. Thus, the pair interaction 
potential in the corresponding boson problem is repul- 
sive, logarithmic in r for r much smaller than A, and falls 
off exponentially as ~ exp(— r/A) when r is larger than 
A. Liquid phases in the vortex system correspond to su- 
perfluid bosonic phases, while vortex crystalline phases 
correspond to a bosonic crystal. 

This mapping between the two systems allows us to 
draw certain conclusions about the behavior of the two- 
dimensional bosonic system from the wealth of informa- 
tion available from existing experimental^— and nu- 
merical^— studies of the mixed phase of type-II lay- 
ered superconductors (high-temperature cuprate super- 
conductors in particular) in the presence of random 



columnar pinning centers oriented parallel to the mag- 
netic field which is perpendicular to the layers. These 
studies establish that when the concentration of pinning 
centers is small compared to that of vortex lines, the 
vortex system exhibits a Bose glass (BoG) phase at low 
temperatures, and a vortex liquid phase at high temper- 
atures. In the vortex-boson mapping, the vortex liquid 
phase corresponds to the superfluid phase of the two- 
dimensional bosonic system. Both experiments^ and 
numerical studies^— show that the BoG phase in the 
vortex system has a polycrystalline structure with grain 
boundaries separatin g cr ystalline domains (see, for ex- 
ample, Fig. 2 of Ref.HH and Fig. 1 of Ref. 111). Exper- 
iments also show the existence of a "vortex nanoliquid" 
phas o 26 i 27 near the boundary between the BoG and vor- 
tex liquid phases (see Fig. 2 of Ref. [13) . From direct vi- 
sualization of the flow of transport current in the system, 
it has been established™ that the "vortex nanoliquid" is 
characterized by the simultaneous presence of solid- and 
liquid-like regions in the system. This phase is distin- 
guished from the homogeneous vortex liquid by the pres- 
ence of interconnected "droplets" of vortex liquid caged 
inside a solid matrix formed by other vortices. The ex- 
periments can not determine the positions of the droplets 
of vortex liquid relative to the grain boundaries that are 
known to exist in the BoG phase. This information is 
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provided by our earlier numerical studie a 28 i 29 of the vor- 
tex system. These studies reproduce all the experimen- 
tally observed features and in addition, show that the 
liquid regions in the "vortex nanoliquid" phase lie along 
the grain boundaries (see, for example, Fig. 7 of Ref. [29]). 
Since the vortex liquid corresponds to the superfluid in 
the bosonic system, these results establish that a zero- 
temperature two-dimensional system of bosons, interact- 
ing via a repulsive pair potential proportional to Kq{t) 
and in the presence of a small concentration of attrac- 
tive pinning centers, exhibits a polycrystalline phase with 
superfluidity along the grain boundaries for a suitable 
choice of system parameters. In the present study, we 
examine whether this behavior persists in the bosonic 
system when the pinning potential is turned off, and cal- 
culate how the rotational inertia of the system is affected 
by the presence of droplets of superfluid lying along the 
grain boundaries. 



B. Methods 

Since our results for the bosonic system are obtained 
from studies of a system of vortices using the methods 
of our earlier workSSr—, we provide here a summary of 
the vortex model we consider and the numerical method 
we use to study its equilibrium behavior. The details of 
both the model and the method of calculation may be 
found in Ref. [H. 

The system studied consisted of a set of vortex lines 
in a highly anisotropic, layered superconductor with the 
magnetic field perpendicular to the layers. The vor- 
tex lines are formed by stacks of "pancake" vortices lo- 
cated on the layers. The Josephson interaction between 
pancake vortices on different layers was neglected. The 
structure and thermodynamic properties of the vortex 
system are determined from numerical minimization of 
a model free energy functional of the Ramakrishnan- 
Yussouff forrr>2ii that expresses the free energy of the sys- 
tem as a functional of the time-averaged local density 
of the pancake vortices. For columnar pins perpendic- 
ular to the layers, the pinning potential is the same on 
all the layers. This implies that the time-averaged den- 
sity distribution is also the same on all the layers. This 
makes the problem effectively two-dimensional 35 . Differ- 
ent phases of the vortex system correspond to different 
local minima of the free energy and phase transitions are 
signalled by crossings of the free energies of different local 
minima. From the density distribution at a local mini- 
mum representing a particular phase of the system, the 
positions of the vortices in that phase are generated by 
locating the points at which the density exhibits local 
peaks. This allows the study the real-space structures of 
different phases. The heights of the local density peaks 
are used to determine whether the vortices in a region of 
the sample are in a solid or liquid state. 

Since our method of calculation is designed for the vor- 
tex system, the parameter that we control is the temper- 



ature T of the vortex system. Changing this temper- 
ature is equivalent to "changing the value of ft" in the 
corresponding zero T bosonic system via the mapping 
Eq. (|2.4I) . This should be interpreted as changing system 
parameters in such a way that the relative importance 
of quantum effects is modified ( "increasing ft" implies in- 
creasing the importance of quantum fluctuations). This 
correspondence should be kept in mind in the interpreta- 
tion of our results in the context of the bosonic problem. 



III. RESULTS 

As mentioned above, all our results for the two- 
dimensional bosonic system have been obtained from nu- 
merical calculations carried out for the equivalent vor- 
tex system. In our earlier studies 2 ^— , we considered 
a system of vortex lines in the presence of randomly 
placed columnar pinning centers and found polycrys- 
talline BoG states when the areal density of the pinning 
centers is much smaller than that of the vortex lines. 
Since quenched disorder arising from the presence of pin- 
ning centers is not present in the 4 He samples studied 
experimentally in the context of supersolid behavior, we 
first investigated whether the polycrystalline states found 
in our earlier studies survive when the random pinning 
potential is turned off. 

We use for this purpose two of the vortex configura- 
tions obtained and studied in Ref. H(J These correspond 
to results for a relative concentration c = 1/8 of columnar 
pins (c = n p /n v where n p , n v are, respectively, the areal 
densities of pinning centers and vortex lines) in a sys- 
tem that contained 4096 vortices (therefore 512 columnar 
pins). We consider, for our starting point, vortex config- 
urations at a vortex temperature T — 17.0K at which 
point (for c = 1/8) the vortex system is well into the 
BoG phase, see the phase diagram in Fig. 1 of Ref. 
(to set the temperature scale, we note that when no pin- 
ning is present, the vortex system considered undergoes a 
first-order melting transition from a crystalline Abrikosov 
lattice to a vortex liquid at T = 18AK .) A perfectly 
crystalline state without grain boundaries is the absolute 
minimum of the free energy in the absence of pinning. 
However if, using the polycrystalline BoG configuration 
as a starting point, one reduces the strength of the pin- 
ning potential to zero in sufficiently small steps while 
running at every step the free-energy minimization rou- 
tine, one ends up, when zero pinning strength is finally 
reached, with a polycrystalline sample, with well-defined 
grain boundaries. This polycrystalline state is of course 
metastable, but we find that as it is warmed up (we use 
steps of 0.2 K), the sample remains in the local poly- 
crystalline free energy minimum (within the finite accu- 
racy of our numerical minimization procedure) until near 
T = 18AK when, as mentioned above, the microcrystals 
melt. Even if these polycrystalline states are not true 
local minima of the free energy, they can be made so by 
the introduction of a few pinning centers that are always 
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present in any physical sample. 
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FIG. 1: (Color online) A metastable poly crystalline vortex 
configuration obtained as explained in the text. The temper- 
ature of the vortex system is T = 18.2 K, slightly below the 
melting temperature of the crystalline solid. Distances are in 
units of ao (see text). The top panel shows the results of a 
Voronoi construction (see text) that brings out the details of 
the poly crystalline structure. In this plot, vortices with 4, 5, 
6 and 7 neighbors are shown as (green) triangles, (red) circles, 
(black) dots and (blue) squares, respectively. Adjacent pairs 
of 5- and 7-coordinated sites correspond to dislocations which 
line up along the grain boundaries that separate the micro- 
crystals. The bottom panel shows the degree of localization 
of the vortices in the same configuration. Here, (red) trian- 
gles and (black) dots respectively represent liquid-like (less 
localized) and solid-like (strongly localized) vortices. This 
plot illustrates the phenomenon of premelting along the grain 
boundaries. 

An example of such polycrystalline configurations, ob- 
tained by warming up the polycrystalline state originally 
obtaind at a vortex temperature T = 17.0K to 18. 2K, is 
shown in Fig. [TJ There we plot the average positions of 
the vortices, defined to be the computational lattice sites 
at which the density has a local maximum. In this and 
subsequent figures, the unit of length is ao, defined by the 
relation T:a^n v = 1 where n v is related to the magnetic 



induction B (which was B = 0.2T in the case shown) by 
n v = -B/<£>o where <E>o is the superconducting flux quan- 
tum. In the top panel, we have shown the results of 
performing a Voronoi construction for the vortex configu- 
ration. The Voronoi cell associated with a lattice point is 
the region of space nearer to that point than to any other 
lattice point and the number of sides of the Voronoi cell 
represents the number of neighbors of the lattice point. 
In the top panel, vortices with six neighbors are shown as 
(black) dots and those with five or seven neighbors are 
shown as (red) circles and (blue) squares respectively. 
A neighboring pair of five- and seven-coordinated vor- 
tices constitutes a dislocation and grain boundaries cor- 
respond to arrays of such dislocations. Microcrystals can 
be seen as ordered regions consisting of sites with six 
neighbors, and these microcrystals are separated by ar- 
rays of dislocations (pairs of sites with five and seven 
neighbors) denoting the grain boundaries. 

In the bottom panel of Fig. [T] we show, for the same 
configuration, the spatial distribution of the value of the 
vortex density at each local density peak, normalized to 
the average vortex density value n v . The value of this di- 
mensionless quantity provides a measure of the degree of 
localization of the corresponding vortex. Using the value 
of this quantity, we can determine whether the vortices 
in a given region are liquid- or solid-like. In our ear- 
lier studies of the vortex system^— (see, for example, 
Sec. IIIC of Ref. El or Sec. HID of Ref. 0), we found 
that vortices in liquid regions correspond to local density 
peaks for which this quantity has values of three or less, 
while those with higher values of this quantity correspond 
to solid regions. In the plot, local density peaks repre- 
senting solid and liquid regions according to this quan- 
titative criterion are denoted as (black) dots and (red) 
triangles, respectively. One can see that all the vortices 
in the ordered regions are solid-like: these regions corre- 
spond to microcrystals. In contrast, the vortices near the 
grain boundaries form liquid channels that separate the 
microcrystals. This plot illustrates the occurrence of a 
kind of "premelting" in the regions near the grain bound- 
aries: these regions melt and become liquid-like (super- 
fluid in the equivalent bosonic system) at temperatures 
lower than that at which the bulk crystal melts (which 
we recall is 18AK). Premelting along grain boundaries 
is well-known^ in classical solids. Our observation of 
this phenomenon in the vortex system is consistent with 
experimental results. At lower temperatures, the liquid 
channels along the grain boundaries are narrower and 
obstructed at the narrowest spots (see Fig. [2] below for 
an example) , forming a chain-like structure of small liq- 
uid droplets. As the temperature is reduced further, the 
liquid regions near the grain boundaries disappear com- 
pletely and are replaced by narrow strips of disordered 
solid. 

These results imply, from the vortex-boson mapping, 
that a two-dimensional polycrystalline system of bosons 
exhibits superfluidity along grain boundaries over a suit- 
able range of parameter values. Using the realistic net- 
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works of grain boundaries obtaind from our numerical 
studies of the vortex system, we can then investigate how 
the NCR! arising from the presence of superfluid channels 
in the bosonic system depends on the system parameters. 
The numerically obtained samples in themselves are too 
small to allow a realistic study of the flows, but a suffi- 
ciently, indeed arbitrarily large, sample can be obtained 
from them through periodic repetitions (periodic bound- 
ary conditions are used in our original numerical studies) 
that produce a tiling of a larger region with the smaller 
samples. One obtains then a larger sample for which the 
short-distance structure is the same, but which has a dif- 
ferent large-scale structure because longer liquid channels 
are produced when different copies of the original sample 
are juxtaposed. Indeed, the length of the largest chan- 
nels after such a juxtaposition is of the order of that of 
the combined sample. It is true that the structure of 
the original small sample is preserved in the small-scale 
structure of the large sample, which has an artificial pe- 
riodicity. However, this is not important: the behavior 
of the moment of inertia of the liquid region, and there- 
fore the NCRI is of course determined essentially by the 
largest-scale channels in the problem. Scaling the system 
by a linear factor of s > 1 in this way introduces liquid 
channels that are larger by a factor of up to s than the 
original ones, while leaving the original small-scale struc- 
tures and the channel widths unchanged. It follows from 
geometrical and similarity considerations that the overall 
moment of inertia of the channels, including the contribu- 
tions from the additional smaller structures, should then 
scale as s . Hence, our conclusions for the moments of 
inertia, which we will normalize to the overall rigid-body 
value, which also scales in the same way, will be indepen- 
dent of this procedure. The samples considered in our 
NCRI calculations are selected portions of 3 x 3 tilings 
of the originally obtained samples. 

The evolution of the network of liquid channels in one 
of these samples with increasing temperature T of the 
vortex system is displayed in Fig. [5] Each panel shows 
a different value of T ranging from 17.8K to 18. AK. In 
these plots, vortices in the liquid regions are shown as 
(blue) squares, and those in the solid regions as (black) 
dots. At the lowest temperature, T = 17.8K, the liq- 
uid regions form largely disconnected, small beads lo- 
cated along the grain boundaries, with larger droplets 
forming at the intersections of two grain boundaries. As 
T is increased, the bead-like liquid structures coalesce 
into channels, so that at about T = 18.2K, connected 
structures of size comparable to that of the whole sam- 
ple (in particular, a ring-like structure that can be seen 
in the upper left corner of the figure) appear. Finally, at 
T = 18.4 K, just below the melting temperature of the 
vortex solid, the liquid regions have fully developed into 
open channels which percolate through the sample. As 
the superfluid in the equivalent bosonic system is allowed, 
in the configurations obtained in the higher T range, to 
freely flow along these channels, whose characteristic size 
is now of the order of the sample size, it is clear that the 



moment of inertia will exhibit a substantial reduction, 
leading to appreciable values of the NCRIF. 

To numerically calculate the moment of inertia of the 
bosonic samples, we proceed as follows: we imagine the 
microcrystalline sample rotating about its center of mass, 
with angular speed f2. The solid regions of the sample, 
composed of the microcrystals, rotate of course as a rigid 
body with the angular velocity fl. The material in the 
liquid (superfluid in the bosonic system) channels will 
flow according to a pattern that can be determined by 
numerically solving the fluid flow equations of an incom- 
pressible, irrotational (super)fluidi 33 ! 34 These equations 
can be solved to obtain the velocity field v(r) in the fluid 
regions. The velocity field at the boundaries of the fluid 
part of the sample must satisfy the boundary condition 33 
that follows from assuming that the fluid is confined by 
the rigid walls. This implies that the component of the 
velocity field v(r) along the outward normal n at any 
point on the boundary must be equal to the component 
of the rigid-body velocity fi X r along n at that point: 

v(r) n = (n x r) n (3.1) 

where r is a vector from the center of rotation to a point 
on the liquid region boundary. In addition, if the fluid re- 
gion is not simply connected, then the circulation of the 
velocity field along a closed path surrounding each inner 
boundary of the region must be specified in order to make 
the problem well-defined 3 ^. We assume that the sample 
is accelerated from rest at zero temperature, so that the 
values of these circulations are zero 3 ^. Once the velocity 
field is obtained, the resulting angular momentum, and 
hence the moment of inertia can be straightforwardly ob- 
tained by numerical integration. 

Analytic solution of this problem is feasibl o 33 i 34 only 
for reasonably simple geometries: for the very irregular 
geometries involving the multiple channels and regions 
under study here solution can only be undertaken nu- 
merically. The numerical method we have used involves 
spatial discretization and iterative relaxation 4 ^ to obtain 
a solution of the Laplace equation (see below) with appro- 
priate boundary conditions. There are certain technical 
difficulties involved in doing this. The main one is that 
our definition of whether a part of the sample is solid or 
liquid is based on the values of the density at the local 
peaks which form a discrete lattice. This lattice must 
therefore be our computational lattice for the purpose of 
calculating the velocity field and the moment of inertia 
and we cannot truly gain additional precision by adding 
more points. Any attempt to do so would necessar- 
ily involve some non-controllable interpolation scheme. 
This unavoidable situation leads to some uncertainties in 
our numerical results. In particular, the final results for 
the rotational inertia include some numerical uncertainty 
arising from this minimum spacing of our computtaional 
lattice. Another technical difficulty is that the bound- 
ary condition for the velocity field, Eq. (|3.1[) . is of the 
Neumann form if one uses the standard method 3 — ^ to 



7 



T=17.8K 




-20 20 60 100 140 180 

x 



T=18.0K 




-20 20 60 100 140 180 

x 



T=18.4K 




FIG. 2: (Color online) Evolution of the network of liquid channels with the temperature T of the vortex system in one of the 
samples studied. Each of the four panels corresponds to a value of T as indicated. Vortices in liquid regions are shown as (blue) 
squares and those in the microcrystal are shown as (black) dots. One can see that as T increases, the liquid channels become 
connected and closed loops of size comparable to that of the sample are formed. Units of length are as in Fig. [1] 



express the incompressible velocity field v(r) as 

v x (r) = -d*(T)/dy, v y (r) = d$(r)/dx, (3.2) 

in terms of a stream function \&(r). The simplication 
of replacing the boundary condition of Eq. (|3.1[) by the 
Dirichlet boundary condition, ^(r) = |fir 2 , at all points 
on the boundary, which was used in earlier studie o 33 ' 34 , 
works for some simple geometries, but leads to incorrect 
results when unobstr acted ring-like channels are present. 
This can be easily checked by applying this method to 
determine the velocity field in a superfluid confined be- 
tween two concentric cylinders rotating about their com- 
mon axis. This forces us to use the original Neumann 
boundary conditions of Eq. ()3.1j) . These are more diffi- 
cult and awkward to implement numerically, as they in- 
volve the outward normal n at the boundary, which can 
be defined in different ways when the boundary consists 
of a set of discrete points. With the possibility of reduc- 
ing the computational lattice spacing not being available 
to us, this introduces additional numerical uncertainty, 



as the final results are somewhat sensitive to the way in 
which n is defined. However, these uncertainties do not 
seriously affect the quantitative conclusions that we can 
draw from our calculation, as we shall see below. We 
recall that using a scalar potential for the velocity field 
leads to other undesirable problems 3 ^ even in some ana- 
lytic cases. 

The irrotational nature of superfluid flow in the ab- 
sence of any superfluid vortex implies that V x v = 
in the superfluid regions. This condition and the defini- 
tion of the stream function ^> in Eq. (|3.2[> imply that ^> 
must satisfy the Laplace equation, V 2, 5(r) = 0, inside 
each superfluid channel. Using a triangular computa- 
tional grid of spacing equal to the average inter-particle 
distance and representing the Laplace operator by sym- 
metric differences^, the Laplace equation can be reduced 
to a set of linear equations satisfied by the values of ^ 
at the computational grid points in the interior of each 
superfluid region. The boundary conditions of Eq. (l3.1[) 
can also be written as a set of linear equations involving 
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the values of \t at the grid points on the boundary of a 
supcrfluid region and their nearest neighbors. This set of 
coupled linear equations is numerically solved using iter- 
ative relaxation^. This yields the values of 'F at the grid 
points, from which the velocity field and the rotational 
inertia can be easily obtained. 

In Fig. |3l we present what is the main result of these 
computations for the NCRI. The quantity plotted there 
(circles) as a function of temperature T (which, we reiter- 
ate, is the temperature at which the calculations for the 
vortex system were performed - it should not be confused 
with the temperature of the equivalent bosonic system 
which is zero) is the NCRIF, defined as the difference 
between the moment of inertia of the entire sample as a 
rigid body and its actual moment of inertia when the su- 
perfluid flow in the liquid portions is taken into account, 
normalized by the total rigid-body moment of inertia. 
These results are for the sample shown in Fig. [3] At 
lower temperatures, when the liquid regions form small, 
disconnected clusters (see the discussion above in con- 
nection with Fig. [5]), the superfluid flow is negligible and 
the entire system behaves as a rigid body, so the value 
of the NCRIF is very close to zero. At higher tempera- 
tures, when the disconnected liquid regions join together 
to form extended channels, superfluid flow begins to oc- 
cur over relatively large portions of the sample and even- 
tually the characteristic channel length becomes of the 
order of the system size. At T = 18.2K, a large, ring-like 
channel opens up (see Fig. [5]) , providing a closed path for 
the flow of the superfluid. This drastically reduces the 
moment of inertia; recall for example that the moment of 
inertia of a superfluid confined in a ring that is rotating 
about its center is zero. Therefore the NCRIF increases 
very drastically (note the logarithmic vertical scale) and 
by the time the temperature is reached where the chan- 
nels are fully open, it has increased by nearly three orders 
of magnitude, while of course remaining relatively small 
(~ 5%). Our conclusion about a rapid increase in the 
value of the NCRIF at T = 18.2K, when a large, ring- 
like channel opens up to allow the flow of the superfluid 
along a closed path, is not affected by the above men- 
tioned uncertainties in the numerical results, which are 
^25% in the worst case. We have repeated the whole cal- 
culation for a second polycrystalline configurations and 
obtained very similar results. Specifically, we find in both 
cases that the NCRIF increases by a factor of ~ 18 at 
T = 18.2K when a large, ring-like liquid channel opens 
up. 

In Fig. |3l we have also shown (triangles) the depen- 
dence of /, the fraction of the sample that is liquid-like 
(superfluid in the bosonic system), on T. It is clear from 
this plot that the dependence of this quantity on T is 
qualitatively different from that of the NCRIF: the liq- 
uid fraction / increases smoothly with increasing T and 
does not exhibit the rapid increase seen in the NCRIF 
at T = \8.2K . This observation brings out the impor- 
tant point that the NCRIF may not provide a good mea- 
sure of the fraction of superfluid regions of the system 



if the superfluidity occurs along a network of channels - 
the connectivity of the channels plays an important role 
in determining the value of the NCRIF. The NCRIF re- 
mains small as long as the network does not contain large 
ring-like paths along which the superfluid can flow with- 
out any blockage. The first opening up of such paths as 
some system parameter is varied leads to a large increase 
in the value of the NCRIF. This large increase, which may 
look like the onset of superfluidity, does not necessarily 
correspond to a sudden, large increase in the superfluid 
fraction / of the system. This should be kept in mind 
while interpreting experimental data for the NCRIF. 
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FIG. 3: The NCRIF, calculated and normalized as explained 
in the text, plotted as a function of T, the temperature of the 
vortex system (circles). Also shown are the values of /, the 
liquid (superfluid in the bosonic system) fraction (triangles). 
The straight line segments join consecutive data points. 



IV. SUMMARY AND DISCUSSIONS 

In this work, we have studied the NCRI arising from 
superfluidity along grain boundaries in a two-dimensional 
bosonic system at zero temperature. By making use of 
standard mappings, we have related the properties of the 
bosonic system to those of a system of superconducting 
vortex lines at finite temperatures. The latter system has 
been extensively studied by numerical methods and com- 
putationally obtained structures in the different phases 
it can exhibit are available. One of the possible struc- 
tures of this vortex system is a polycrystalline solid in 
which liquid regions near the grain boundaries between 
microcrystals expand and form more connected regions 
which eventually extend through the sample as the vortex 
system temperature is increased. These vortex configu- 
rations, obtained from previous numerical studies, cor- 
respond, via the above mentioned mappings, to different 
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configurations of the boson system, each containing more 
or less extensively connected narrow superfluid regions, 
separating much larger and compact crystalline regions. 
We have then calculated the moment of inertia of the 
bosonic samples by numerically solving the equations of 
superfluid hydrodynamics in the superfluid regions of the 
sample. 

We have found that the non-classical portion of the 
moment of inertia of such samples (the NCRIF) is nearly 
zero as long as the vortex liquid (bosonic superfluid) re- 
gions remain disconnected, but it increases abruptly by 
about two orders of magnitude as the superfluid chan- 
nels become connected to form ring-like structures with 
sizes comparable to the system size. The values we find 
for the fully developed NCRIF are of the order of a few 
percent. Thus, our results indicate that NCRI behav- 
ior can indeed arise in two-dimensional bosonic systems 
from superfluid regions associated with grain boundaries. 
The magnitude of the jump in the NCRIF and the final 
NCRIF values obtained from our calculations are quite 
similar to what is found experimentally in solid 4 He. Our 
results, therefore, lend support to the notion that super- 
solid phenomena in 4 He are related to superfluidity in 
grain boundaries or other crystal defect regions, such as 
dislocations. Our results also indicate that an abrupt 
increase in the NCRIF from a vanishingly small value 
to a few percent does not necessarily correspond to a 
similar increase in the fraction of the sample that is su- 
perfluid. As shown in Fig. [3J such an increase in the 
NCRIF may actually correspond to the opening up of 
system-spanning ring-like channels through which the su- 
perfluid can flow without any block. The importance of 
the existence of a percolating network of superfluid chan- 
nels in observing macroscopic signatures of superfluidity 
has been pointed out in an earlier study^i. However, we 
are not aware of any calculation of the behavior of the 
NCRIF across the percolation transition. 

As discussed above, the results of existing experimen- 
tal^r— and numerical^ - — studies of vortex lines with 
dilute columnar pinning strongly suggest, via the map- 
ing used here, that a polycrystalline state with superflu- 
idity along grain boundaries occurs in a two-dimensional 
bosonic system with a small concentration of strong pin- 
ning centers. Such a state may be seen in experments 41 
on films of 4 He adsorbed on substrates with imperfec- 
tions that act as pinning centers. In the vortex system, 
the equilibrium phase changes from polycrystalline Bose 
glass to vortex nanoliquid and then to homogeneous liq- 
uid as the temperature is increased. In the equivalent 
bosonic system, this sequence corresponds to going from 
a polycrystalline solid to a "supersolid" with superfluid- 
ity in the grain boundary regions and then to a homo- 
geneous superfluid. It is, however, not clear how the ex- 
perimentally accessible parameters in the bosonic system 
can be changed to simulate the effects of increasing the 
temperature in the vortex system. The same sequence 
of phases can also be seer*2£r— in the vortex system by 
increasing the magnetic induction at a fixed, low temper- 



ature. Since the magnetic induction in the vortex system 
determines the areal density of the vortex lines, the su- 
persolid phase in a two-dimensional bosonic system may 
be accessed by changing the coverage of the 4 He film. 
The supersolid would be the true equilibrium phase of the 
system when random defects are present. Our numerical 
results, furthermore, suggest that the supersolid survives 
as a metastable phase even in the absence of random 
pinning. While it is not obvious how such a metastable 
phase could be accessed in experiments, the supersolid 
behavior observed in three-dimensional solid 4 He must 
also be a metastable phenomenon if it arises from super- 
fluidity in a network of structural defects, because these 
defects would not be present in the true equilibrium solid 
which should be a perfect crystal. 

In making any comparison of the behavior of the 
bosonic system found in our study with experiments on 
He samples, one should keep in mind the important fact 
that the interaction between two He atoms (strongly re- 
pulsive at short distances, with a weak attractive part at 
larger distances that falls off with distance as a power- 
law) is substantially different from that between two 
bosonic particles in the system considered here (purely 
repulsive, logarithmic at short distances and exponen- 
tially decaying at longer distances). For this reason, the 
results of our study can not be applied quantitatively to 
experimentally studied 4 He systems. However, we expect 
the qualitative behavior found in our study to be ob- 
served in systems of 4 He atoms because the phenomena 
on which our conclusions are based (i.e. the formation of 
polycrystalline structures and "pre-melting" along grain 
boundaries) are fairly generic, independent of the details 
of the inter-particle potential. 

In experiments on 4 He, a normal solid to supersolid 
transition is observed— ~— as the temperature of the sam- 
ple is decreased. So, it is interesting to inquire whether 
our study provides any information about the behavior of 
the bosonic system as its temperature is changed. Unfor- 
tunately, we can not address this experimentally relevant 
question in our study. This is because a nonzero temper- 
ature in the bosonic problem maps to a vortex system of 
finite thickness L, and our studies of the vortex system, 
carried out for the L — > oo limit, do not provide any in- 
formation about phase changes that may occur as L is 
varied at constant temperature. 

We close with a discussion of whether the available ex- 
perimental results for supersolid 4 He show any evidence 
for the sequence of phases (polycrystal to vortex nano- 
liquid to vortex liquid) found in the vortex system upon 
increasing its temperature T. In the equivalent bosonic 
system, this sequence correspond to a transition from a 
defected solid to a supersolid, followed by a second tran- 
sition from the supersolid to a superfluid. The vortex- 
boson mapping implies that increasing T in the vortex 
problem is equivalent to increasing the relative impor- 
tance of quantum effects in the zero-temperature boson 
problem. It has been suggeste d 42 i 43 that increasing the 
pressure P has the effect of reducing the relative impor- 
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tance of quantum fluctuations in solid 4 He. Assuming 
this to be correct, increasing (decreasing) T in the vor- 
tex system should be analogous to decreasing {increas- 
ing) P in 4 He experiments. Then the relevant question 
is whether a sequence of normal solid to supersolid to 
superfluid transitions occurs in 4 He as P is decreased at 
constant temperature. There is some experimental evi- 
dence 4 ^ for a reduction in the low-temperature value of 
of the apparent superfluid fraction in supersolid 4 He as 
P is increased beyond 55 bar. This may correspond to a 
transition from the supersolid to a regular defected solid, 
analogous to the disappearance of the liquid regions in 
the vortex system as T is decreased. Also, the "phase 
diagram" of 4 Hc in the P — T plane shown in Ref. G] 
suggests that the temperature of the normal solid to su- 
persolid transition decreases as P is increased (the tran- 
sition temperature is found to decrease from 315 mK at 



a pressure of 26 bars to 230 mK at pressures exceeding 
40 bars) . If this relatively weak dependence of the tran- 
sition temperature on pressure is a genuine effect, then 
decreasing P at constant T would indeed lead to a se- 
quence of two transitions, first from the normal solid to 
the supersolid and then from the supersolid to the su- 
perfluid phase. This would be analogous to the behavior 
seen in our calculations. More detailed investigations of 
the pressure-dependence of supersolid behavior in He 
would be very interesting in this context. 
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